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Based on calculations of the electronic structure of coupled multiple quantum dots, we study 
systemically the transport properties of the system driven by an ac electric field. We find quali- 
tative difference between transport properties of double coupled quantum dots (DQDs) and triple 
quantum dots. For both symmetrical and asymmetrical configurations of coupled DQDs, the field 
can induce the photon-assisted Fano resonances in current-AC frequency curve in parallel DQDs, 
and a symmetric resonance in serial DQDs. For serially coupled triple quantum dots(STQDs), it 
is found that the A-type energy level has remarkable impact on the transport properties. For an 
asymmetric (between left and right dots) configuration, there is a symmetric peak due to resonant 
photon induced mixing between left /right dot and middle dot. In the symmetric configuration, a 
Fano asymmetric line shape appears with the help of "trapping dark state". Here the interesting 
coherent trapping phenomena, which usual appear in quantum optics, play an essential role in quan- 
tum electronic transport. We provide a clear physics picture for the Fano resonance and convenient 
ways to tune the Fano effects. 
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I. INTRODUCTION 



The electronic transport in quantum dot structures is 
of great interest and has been the subject of active re- 
search in recent years. The strong spacial confinement 
leads to the discrete energy spectrum. These quantum 
dots (so called "artificial atoms" ) /mutiple coupled quan- 
tum dots (MCQDs, so called "artificial molecules" ) are 
very important in understanding fundamental quantum 
phenomena. Many interesting phenomena which appear 
in real atoms and molecules can be well manifested in 
these nanostructurcs. Tunability provides much more op- 
portunities for exploring richer physics which is hard to 
access in real atoms/molecules. It also leads to various 
applications. For instance, novel quantum logical gates 
and elementary qubits in quantum computersii^i^ could 
be realized based on these systems. 

In these nanostructurcs, quantum coherence and in- 
terference is one of the most important issues and leads 
to a lot of interesting phenomena such as Aharonov- 
Bohm (AB) oscillations^, Kondo effect^'^, coherent 
trapping^'®, Fano resonance^' and so on. In re- 
cent years, double and triple quantum dots have at- 
tracted much attention due to their rich electric struc- 
ture and physical phenomena. An AB interferometer 
containing two QDs has been realize d^^'^'^i^'^'^^ . Some 
experimental groups have been able to fabricate triple 
quantum dots with hight qualit y and many rele- 
vant theoretical studies have been conducted on these 
systemsiaa9>20^^^^^^^^^. In spite of many 
studies on the multiple quantum dots, relatively little at- 
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Serial Triple Qantum Dots 

FIG. 1: Schematic diagrams of coupled multiple quantum 
dots 



tention has been paid to the photon-assistant transport 
in multiple quantum dots driven by a time-period electric 
field. Moreover, periodic driven triple quantum dots with 
the relevant three-level structure may have some interest- 
ing interference phenomena, such as coherent trapping, 
like those in quantum optics. In this paper we study the 
photon-assistant transport in MCQDs (see Fig.l) paying 
special attention to the consequences of those interesting 
interference phenomena. 

In previous theoretical work on MCQDs, the quantum 
levels, which depend on the structures of the system, are 
often assumed as parameters^2i3i^ and the quantum prop- 
erties of MCQDs can not be presented quantitatively. 
Therefore, it is necessary to reveal the quantum behav- 
ior of MCQDs based on more realistic model. In our ap- 
proach, we first design DQDs and STQDs (with A-type 
three-level structure) using a two-dimensional confining 
model in the effective mass frame. Based on the obtained 
level structure and with the help of Floquet theory^, we 
study the transport properties of the DQDS and STQDs. 
It is found for both symmetrical and asymmetrical config- 
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FIG. 2: (a) Expected A-type three-level scheme for STQDs. 
(b) Confining potential model of STQDs. (c)The three solid 
lines are levels for a STQD as functions of hujr for STQDs 
of hun — 3.00 meV, hum ~ 3.52 meV, and d = 77.70 nm. 
The dotted two lines are levels for DQDs of hui =3.00 meV, 
d = 77.70 nm. Insert: The eigenwavefunctions for STQDs 
(up) and DQDs (bottom) at hur = 3.00 meV and 2.70 meV 
respectively. 



urations of coupled DQDs, there is photon-assisted Fano 
resonances in parallel DQDs, and a symmetric resonance 
in serial DQDs. The electronic transport is very differ- 
ent from that in STQDs with A-type energy level: the 
photon-assistant tunneling leads to the symmetric Breit- 
Wigner— resonance when the system is asymmetric (be- 
tween left and right quantum dots). When the system 
is in a symmetric (between left and right quantum dots) 
configuration, the formation of "trapping dark state" re- 
sults in the interesting asymmetric Fano resonance under 
resonant condition (photon energy equals to the energy 
difference between left /right quantum dot and the mid- 
dle dot.) Our studies show that the transport properties 
are quite sensitive to the number of quantum dots in the 
coupled systems. 

The organization of the paper is as follows. In sec- 
tion II we describe the model and calculate the electronic 
structure of the MCQDs. In section III we present the 
photon-assistant transport properties of the system and 
discuss the results. A brief summary is given at the end 
of the paper. 



II. MODEL AND ELECTRONIC STRUCTURE 

We use a two-dimensional model to describe multiple 
coupled quantum dots laterally. Here, we give serially 



coupled triple quantum dots as an example. The confin- 
ing potential is 



V{x, y) — 0.5m*min{cLif^(a; + (£f + ^lyV^ , ^mx^' 

I 2 2 2 / i\2 I 2 21 



(1) 



where d is the interdot distance, m* is the effective mass 



of electron. 



and 



are confining trap 



frequencies of the left, middle and right dots in the x{y) 
direction respectively. The model Hamiltonian of an elec- 
tron in the coupled quantum dots can be written as 



H 



2ra* 



V{v), 



(2) 



where r = {x,y). In our calculation, the material pa- 
rameter of GaAs QDs is used as m* = 0.067TOe and the 
value of d is taken to resemble experimental systems. We 
use the eigenstates (^i, , (pi^ , (pi^ for each dot as the ba- 
sis of the Hilbert space. Considering the nonorthogo- 
nality of the basis states, we obtain the eigenstates of 
Eq. (2) by solving the generalized eigenvalue of the sys- 
tem. Here, we investigate the levels of the STQDs by 
varying hurx — hujry — hiUr with constant parameters 
of hwix = hujiy = h^i — 3.00 meV, hujmx = fv-^my = 
hiOm = 3.52 meV, and d = 77.70 nm. For such param- 
eters, the tunneling energy between the left and middle 
dot is about 51 /leV. The lowest three energy levels of 
STQDs form a A-type structure as shown in Fig. 2(a) 
(c). Experimentally, the levels of each dot are tuned 
by changing the voltage of the corresponding electrode 
gate. Such change of voltage corresponds to the vari- 
ation of confining strength of the dot in our model. As 
the value of huJr increases, the right level increases. When 
hwr = hjji = 3.00 meV, the energies of |1) and |3) anti- 
cross. The corresponding eigenwavefunctions are shown 
in the insert. It is clear that a pair of delocalized bond- 
ing and antibonding states are formed, while the other 
state is localized in the middle quantum dot. When hwr 
is away from 3.00 meV, the three levels are nearly en- 
ergies of the ground states of the left, middle and right 
dots respectively, and the eigenstates are all localized in 
each quantum dots. The continuous manifolds on the 
two sides of the Fig. 2(a) correspond to electronic states 
with chemical potential and and the two side dots 
are coupled to the leads with the dots-lead hopping rate 
Vl and Vr- 

Using the same method, we investigate the levels of 
the DQDs by varying hwr with constant parameters of 
fvjji = 3.00 meV, and d — 77.70 nm. The lowest two 
energy levels and the corresponding eigenwavefunctions 
with different values of hwr of DQDs as shown in Fig. 2(c). 
Similarly, when hwr — tiwi = 3.00 meV, the energies of 
|1) and 1 2) anticross and a pair of delocalized bonding 
and antibonding states are formed. When tujjr is away 
from 3.00 meV, the two levels are nearly energies of the 
ground states of the left and right dots respectively, and 
the eigenstates are localized in each quantum dots. 
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III. PHOTON-ASSISTANT TRANSPORT 
PROPERTIES 

A. Formulism for photon-assistant transport in 
coupled QDs system 

We first construct our formulism based on parallelly 
coupled DQDs (see Fig.l ) with an external time- varying 
field. This system can be described by the following 
Hamiltonian: 



H= Ha + Ho + Hj 

a=L,R 



with 



Ha = ^ ^a.k^'a.k^' 



Hd^Y E,(t)d\d, ~ {t,d\d2 + H.c. 

1=1,2 



H' 



T= ^ taid\aa,k + H.C. 



(3) 



(4) 



(5) 



(6) 



Ha {a = L,R) describes the left and right normal metal 
leads. Hd models the parallel-coupled DQD where d] 
(di) represents the creation (annihilation) operator of the 
electron with energy Ei in the dot i {i = 1,2). tc denotes 
the interdot coupling strength, which can be obtained 
from the calculations of the electronic structure. Under 
the adiabatic approximation, the external ac electric field 
can be reflected in the single-electron energies which can 
be separated into two parts as Ei it) = i?i + Ao [t) and 
E2{t) — E2 — Ao(<) for the central conductor. Ei is the 
time-independent single-electron energies without the ac 
field, and Ao(t) is a time-dependent part from the ac 
field, which can be written as Ao(i) = edA cos fit. Ht 
represents the tunneling coupling between the DQD and 
leads where tai is the hopping strength between the ith 
QD and the a lead. To capture the essential physics 
of photon-assisted Fano resonance, we consider the sim- 
plest case with noninteracting electrons in two single- 
level QDs. 

The current la (t) from the a lead to the central region 
can be calculated from standard NGF techniques, and 
can be expressed in terms of the dot's Green function 



as 



33.34.35 



Ia{t) 



2e 
li 



G<{t,t') = i(*t(t')\]/(i))^ respectively, with the opera- 
tor ^-t = (dj,4)- 

The total self-energy is S** = SJ^ + SJJ, in which 
SJ^ and SJJ are caused by the a lead and the interdot 
couplings, respectively. Under the wide-band approxi- 
mation, the retarded self-energy caused by the a lead is 
defined as 



^a{t,t') 



-m-t') 



\ /"pQ"pa 

V V 1^ 2 



pa 
^ 2 



(8) 



where Ff is the linewidth function defined by F" = 
2'Kpat*^itai with pa being the density of states of the cor- 
responding lead, describing the coupling between the ith 
QD and the a lead. The advanced self-energy can be 
obtained from '^'i{t,t') = (S^(t,t'))t. The self-energy 
caused by the interdot coupling is 







-tr 



-tc 



(9) 



The lesser self-energy is S<(ii,i2) = X^a ^q(*1j*2) and 



e(t-t') 



11 V^1^2 )Q0) 



where Ja{f) — 1/(6*^ '^"^"^^'^^^ -I- 1) denotes the Fermi 
distribution function of electrons in the a lead. 

The Green's function of uncoupled DQD without the 
couplings to the two leads can be easily obtained as 

After taking the Fourier transformations, the Dyson 
equation and Keldysh equation become 

GL„W = g™„(e) + E GLi(e)S[;(e)gL(e), (12) 



and 



:G;„,(6)S<(6)Gr„(e), 



(13) 



J 

Re J dt'TT{[G^it,t')i:<{t',t)+G<{t,t')i:^^it\t)]}, 

(7) ' 

Here, the Green's function C'^ and the self-energy 
S"'^ are all two-dimensional matrices for the DQD sys- 
tem. The bold-faced letters are used to denote ma- 
trices. The retarded and lesser Green functions are 
defined as G''{t,t') = -iOit - t'){{^{t),^^t')}) and 



respectively. With these Green's functions, the time- 
dependent current can be expressed as 

lait) ^ ^ReEe'""{ / ATr[G'(6)S<(6)+G<(6)S^(6)],o}. 
I 

(14) 

Then, the average current is 



2e, 



Re{ / deTr[G'-(6)S<(e)+G<(e)SS(e)]oo}. 

(15) 



Here we give another formulism suitable for the seri- 
ally coupled QDs system. Making use of the Floquet 
decomposition, the physical picture is clearer in this for- 
mulism. As seen in Fig. 1, we should have the serially 
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FIG. 3: (a) Average current / of the serial DQDs. (b) Average 
current I of the parallel DQDs where hwr = 3.10 meV, fcsT = 
0. (c) (d) are the same as (a) (b) respectively except that 
huor ~ 3.00 meV. 



coupled DQD, if we set r}j = Tj^ — 0. In this case, the 
self-energies are simplified. It can be shown that 

+ 00 „ 

E j de{T'il{s)fn.{e)~T'il{e)fde)}, (16) 



k— — oo 



where 



rl^)(s) = rir;,|(i|GW(e)|iV)p (i7) 



T'il{e)=TnVLmG^'\e)\l)\' (18) 

(N=2) denote the transmission probabilities for electrons 
from the right lead and from the left lead, respectively. 
fi{e) = (1 + exp[(e — iJ,()/kBT])~^ denotes the Fermi 
function and 



|M;3,fc' + fc)(M^,fc/| 

(e/3 + k'Ml - ijfi) 



(19) 



is the Fourier coefficients of the retarded Green function, 
where |u/3,fc) are the Fourier coefficients of the Floquet 
state \up{t)), tp, {(3 = 1, 2) are the real and imaginary 
parts of quasi energies respectively. This formulism is 
valid for general N serially coulped QDs systemi^^ 



B. Transport properties of DQDs 

In our numerical calculations, we assume ksT — 
and set the energy independent dots-lead hopping rate 
Tl — Tji = 9 /ieV, the applied voltage /i^ — fiR — 15 /iV, 
and the ac field magnitude A = 1.54 V/cm. First we 
study the transport properties of an asymmetric system 
by applying different confining potentials to the left and 
right dots. The average current-frequency curves of seri- 
ally and parallel DQDs with hur = 3.10 meV are shown 



in Fig. 3 (a) and (b) respectively. We find that the cur- 
rent curve has a symmetric Breit-Wigner line shape when 
the two dots are serial, and the current has an asym- 
metric Fano line shape when the dots are parallel. Here 
the Breit-Wigner line shape is a consequence of photon- 
assistant resonant, i.e., hft — y/ {Ei — £2)"^ -\- 4t^. The 
new photon-assistant Fano lineshape comes from the 
interference between photon-modified bound-antibound 
channels. Then we apply the same confining potentials 
to the two dots and research the electron transport of a 
symmetric system. The results with Hujr = 3.00 meV are 
displayed in Fig. 3 (c) and (d). Similarly, we find that 
a Breit-Wigner resonance appears in the serially DQDs 
and a Fano resonance occurs in the parallel DQDs. From 
above discussion, we conclude that the ac electric field 
can induce the photon-assisted Fano resonances for both 
symmetrical and asymmetrical parallel configurations of 
DQDs, but can not induce Fano resonance in the serially 
DQDs, whether the system is symmetric or not, which 
means that the energy level of DQDs does not have re- 
markable impact on the transport properties. 



C. Transport properties of STQDs 

As one observes in the last section that there is no 
Fano effect in a serially coupled DQD due the lack of 
apparent interference channels. One may expect that 
there is also no Fano resonance in a STQD. However, our 
studies shows the transport in STQD is quite different 
from that in a DQD, and we will give a quite different 
story in this section. 

Based on the results of electronic structure in section 
II, we study the transport properties through the sys- 
tem with A-type three-level structure ( Fig. 2 ) under the 
action of an ac driving field with a frequency fJ. 

In our numerical calculations, we set the applied volt- 
age /zl — /i_R = 120 /LtV, and the other parameters we have 
chosen are the same as the ones for DQDs. 



We first study the transport properties of an asymmet- 
ric system by applying different confining potentials to 
the left and right dots. The average current / as a func- 
tion of the driving frequency U, with fiWr = 3.10 meV 
is presented in Fig. 4 (a). We find the current curve 
has a symmetric Breit-Wigner line shape around il — 
129.48 GHz, suggesting that there is a resonance for 
electrons in the system in this case {h£l — E2 — Ei). 
The time evolution of the probabilities for an electron 
in left, middle, right dot are shown in Fig. 4 (b). Here 
we have performed our calculation in a closed system 
(i.e., without interaction with the leads) and have used 
the initial condition = 1 and the resonant condi- 
tion Vl = 129.48 GHz. Fig. 4(c) shows the time av- 
erage of the probability distribution of Floquet states, 
Pp = ^ \up{t)\'^dt, /3 = 1, 2, 3. It is clear that there is 
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FIG. 4: (a) Average current J as a function of the driving frequency 57. (b) Time-dependent occupation probabilities for electron 
in left, middle, right dots, (c) Time average of electron probability distributions of the Floquet states sd, Q — 129.48 GHz where 
hur ~ 3.10 meV, fcsT = 0. (d) (e) (f) are the same as (a) (b) (c) respectively except that Q, — 130.57 GHz and hujr = 3.00 meV. 



a photon assistant mixing between left and middle dot. 
This photon assistant charge transfer leads to the occur- 
rence of the current resonance phenomenon. 

Then we consider the case that the system is symmetric 
by applying the same confining potential to the two side 
dots. The average current / as a function of the driving 
frequency with hiUr = 3.00 meV is displayed in Fig. 4 
(d). We find that the current curve has an asymmetric 
Fano line shape around 130.57 GHz and the current 
amplitude is much larger than that in the asymmetric 
case. In order to understand the intriguing phenomenon, 
we did similar calculation of the time evolution of the 
occupation probabilities for an electron in left, middle, 
right dot and the time average of probability distribu- 
tion of Floquet states. As shown in Fig. 4 (e) and (f), 
we find that two delocalized Floquet states ( which is 
related to bonding and antibonding states ) are formed 
and another Floquet state is still localized in the middle 
dot in this case, indicating that the middle dot medi- 
ates the super-exchange interaction between the left and 
right dots. The states |1), |3) behave like "trapping dark 
state" in atomic system. Though state |2) has a very 
small occupation probability in the left and right dots, 
it plays an important role in the occurrence of the Fano- 
type resonance, which can be presented in the following 
discussion. 

As seen from Eq. (19), three Floquet states form three 
conducting channels. Due to the interaction with leads, 
the three states may have different width. In the sym- 
metric situation, the bonding and antibonding Floquet 
states I Ml), Iws) have large occupation probability in left 
and right dots, thus are strongly coupled to the leads. 
They have wide width and can be viewed as continuous 
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FIG. 5: Phases as functions of the driving frequency fl. 



channels. The localized Floquet state \u2) has small oc- 
cupation probability in left and right dots and is weakly 
coupled to the leads. It has very narrow width and can 
be viewed as resonant discrete channel. The interference 
between strongly coupled channels and weakly coupled 
channel leads to the observed Fano resonance22. While 
in the asymmetric situation, the Floquet state |m2) is 
also effectively coupled to the leads via a photon-assistant 
mixing between left and middle dot. Thus all three states 
form strongly coupled channels, so no Fano effect is ob- 
served. Here we discuss in detail the interference between 
different channels. From Eq. (19), one can see that the 
average current / is mainly decided by K-ZVIG'*^^ (e)|l)p. 
The interference effects can be seen from following calcu- 
lation (here we show term with fc = 0, e = 60 as an 
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example) 

(20) 

In Fig. 5, we show the phases 0/3(/3=i.2,3) as the func- 
tions of the driving frequency Vl. Here, the phases 4>i 
and (j)^ are phases of two strongly coupled states and 
M3}, and 02 is the phase of the weakly coupled state \u2)- 
From Fig. 4, we find that has a phase shift that varies 
from to TT as the field frequency f2 is moved through the 
resonance frequency, which provides resonant channel for 
electrons to transport the system, , (^3 arc nearly con- 
stant and can be assumed to be independent of the driv- 
ing frequency, which can be considered as background 
or nonresonant channels. Therefore, the interference be- 
tween the strongly coupled channels and weakly coupled 
channel results in the Fano-type resonance. Finally, we 
would like to point out that the above picture also ex- 
plains why there is no Fano resonance in a serially cou- 
pled double dot. As a matter of fact, in that situation, 
there are only strongly coupled bonding and antibonding 
states. Therefore, there is no apparent Fano effect. 

In summary, we have studied the transport properties 
of DQDs and STQDs under the action of an ac electric 
field. Using a two-dimensional confining potential model 
in the effective mass frame, the two-level structure and 
A-type three-level structure are obtained through solv- 
ing the generalized eigenvalue of the system. Based on 
the level structure and Floquet theory, we investigate the 



dependence of the ac current on the electronic structure 
of the system. It is found that the two-level structure 
does not influence the transport properties of DQDs re- 
markably. For both symmetric and asymmetric configu- 
rations, there is Fano resonance in parallel DQDs and no 
Fano resonance in serial DQDs. However, it is a different 
case for STQDs. The A-type three-level structure has 
great impact on electron transport: when the system is 
asymmetric, the symmetric Brcit-Wigncr resonance ap- 
pears due to phonon assistant tunneling; When the sys- 
tem is symmetric, the interesting asymmetric Fano res- 
onance occurs under resonant condition. In this case, 
quantum interference results in the formation of trapping 
dark states. These trapping dark states are the delocal- 
ized bonding and antibonding states which serve as the 
continuous channels. Here the middle dot plays a dual 
role: (1) It mediates the super-exchange interaction be- 
tween the left and the right dots; (2) The localized state 
within it serves as the resonant discrete channel. Our 
work presents the unique quantum interference features 
of the multiple quantum dots system and is useful for the 
design of novel nanodevices. In short, for transport in a 
system with N quantum dots, the situation with N=l, is 
different from that with N=2 and the case with N=2 is 
different from that with N=3. We do not expect quali- 
tative difi^erence for system with N > 3. 

This work is supported in part by the National Natural 
Science of China under No. 10574017, 10774016 and a 
grant of the China Academy of Engineering and Physics. 
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